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ABSTRACT 

Imaging at 1.25/im and 2.20/xm has been obtained of the field containing 
the galaxy (RDl) found at redshift z — 5.34 by Dey, et al.(1998). This galaxy 
has been detected at 1.25//m, while the lower redshift {z — 4.02) galaxy also 
found in the same field by Dey, et al. was detected at both 1.25/xm and 2.20//m. 
Comparison to stellar population synthesis models indicates that if RDl is a 
young (t < 10^ yrs) galaxy, significant reddening (Ay > 0.5 mag) is indicated. 
Combined with observations of other high redshift systems, these data show 
that dust is likely to be an important component of young galaxies even at 
redshifts of z > 5. The extinction-corrected monochromatic luminosity of RDl 
at 1500A is then a factor of about three larger than L^gQo as determined by 
Dickinson (1998) for 2; ~ 3 starburst galaxies. The implied star formation rate 
in RDl, corrected for extinction, is ~ 50 - 100/150^ Mgyr-^ 

Subject headings: galaxies: individual, evolution; dust, extinction 
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1. Introduction 

In the last several years, the combination of new wavelength dropout discovery 
techniques (e.g. Steidel et al. 1996) coupled with the deep imaging power of the Hubble 
Space Telescope and the spectroscopic capabilities of a new generation of large ground-based 
telescopes, has lead to a blossoming of the study of galaxies at redshifts of ^ > 2. A recent 
breakthrough in this field has been the discovery of a galaxy at a redshift of 2; = 5.34 by 
Dey et al. (1998). Besides being the highest redshift object presently known, this galaxy is 
important because it is apparently neither gravitationally lensed, nor the host of a powerful 
active nucleus. Thus there is hope that detailed studies of this system and others like it at 
2; ~ 5 will shed light on the intrinsic, stellar properties of very young galaxies. 

Because visible observations of objects at such high redshifts are sampling the restframe 
ultraviolet part of the spectrum, they are very sensitive to the presence of even small 
amounts of intervening dust. Dust is important at these redshifts for three simple reasons. 
First, intrinsic luminosities and star formation rates are always underestimated when dust 
is not taken into account. This may have vast implications for metal production at early 
epochs (Madau et al. 1996, Meurer et al. 1997). Second, dust may hide entire populations 
of high redshift, star-forming galaxies from detection because of sample selection effects (e.g. 
Dey, Spinrad & Dickinson 1995). Third, dusty interstellar media in galaxies at 2; > 3 — 4, 
by their very nature, imply previous episodes of star formation and enrichment, pushing 
the first epoch of star formation to even higher redshifts. Although "pollution" of the ISM 
can happen relatively quickly from a truly primordial state (e.g. Matteucci & Padovani 
1993), dusty systems at high redshift always provide a lower hmit to the redshift of the 
earliest star formation. Since there is mounting evidence for dust in high redshift galaxies 
(e.g. Dunlop et al. 1994, van Ojik et al. 1994, Ivison 1995, Dey, Spinrad & Dickinson 1995, 
Knopp & Chambers 1997, Sawicki & Yee 1998) it is natural to ask whether dust is present 



-4- 



in galaxies ai z > 5. 

To understand the true nature of the z — 5.34 galaxy discovered by Dey, et al. (1998), 
hereafter referred to as RDl - for Red Dropout 1, it is necessary to know whether it is being 
viewed through a significant amount of dust. To achieve this, it is essential that such objects 
be observed at the longest possible wavelengths. We report in this paper serendipitous 
observations made of this galaxy at J (1.25//m) and K (2.20 yum). The galaxy is clearly 
detected at J, and simple fits to the observed optical-infrared colors suggest a reddened, 
young stellar population. In the discussion we adopt the cosmological parameters Hq = 50 
km s~^ Mpc~^ and qo =0.5. Throughout this paper we refer to broad-band magnitudes 
measured with respect to Vega, such that J= mag and K= mag correspond to 1578 Jy 
and 646 Jy, respectively. 



2. Observations and Data Reduction 

As part of an observational study of high redshift radio galaxies (Armus, et al. 1998a, 
1998b) we observed the most distant known radio galaxy, 6C0140-I-32 in August 1997 and 
January 1998 using the near infrared camera (NIRC) on the Kl telescope of the W.M. Keck 
Observatory. The instrument is described in detail by Matthews and Soifer (1994). It has 
a 256x256 InSb array with 0.15"x0.15" pixels for a 38"x 38" field of view. On 18 August 
1997, the field was observed at K under conditions of thin cirrus and good seeing (~ 0.5"). 
On January 16, 1998 we observed this field under clear, photometric conditions at J and K 
in excellent seeing (0.4 — 0.5" FWHM). On 17 Jan 98 we also observed the field at K under 
conditions of good seeing with thin cirrus. 

The target for the observations was centered in the field of view, and observations 
were made with individual integration times of 60 seconds. In the August observations the 
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target was moved in a regular 9 point square grid covering a size 10" on a side, with the 
center of the grid moved by 1 — 2" between grids. In the January observations, the target 
was moved randomly over an area of 10" x 10" between individual observations in order to 
affect optimal sky subtraction. 

Because the galaxy RDl and the companion galaxy BD3 (for Blue Dropout 3) at 
a redshift z — 4.022 were only ~ 12" from the targeted object, they were contained in 
virtually every frame, and the combined data achieved the full sensitivity of the observation 
at the location of RDl and BD3 at each wavelength. 

A bright star in the center of the NIRC field was used to accurately register the 
successive frames for coaddition. To remove time-variable fluctuations in illumination, 
separate sky and normalized flat-field frames were created from the data for each set of 
three images, by taking the median of the nearest 7-9 frames. After being trimmed to a size 
of 251 X 251 pixels, the individual data frames were sky subtracted, flat flelded, and shifted 
to a common dc level after known bad pixels were flagged. These processed images were 
then aligned, using integer pixel shifts, and combined using a clipped mean algorithm. The 
FWHM of a point source in the flnal mosaic, as determined from several stellar images in 
the field, was approximately 0.5" at J and 0.45" at K. The photometry was based soley on 
the data from 16 Jan 98, when conditions were photometric. The data were calibrated by 
reference to the infrared standard stars of Persson et al. (1998). 

3. Results 

The J and K images of the field centered on RDl are shown in Figure 1. The J data 

are from 16 Jan 98 only, while the K data are from 18 Aug 97, 16 Jan and 17 Jan 98. We 
obtained a total of 70 and 92 minutes of integration at J and K, respectively. 
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In order to bring out the faint, high redshift sources RDl and BD3, the data in Figure 
1 have been convolved with a circular gaussian having a FWHM equal to the measured 
stellar FWHM in the final J and K-band mosaics. In this figure, 12" sections of the mosaics 
are displayed. Objects RDl and BD3 are both detected in the J-band image. At K, BD3 is 
marginally detected, while RDl is undetected. The J-band image of RDl is resolved in the 
north-south direction, while that of BD3 is unresolved. 

Because the objects RDl and BD 3 are separated by ~ 1.2" obtaining photometry of 
both systems required special care. A square beam, 0.75" on a side, was centered on both 
RDl and BD3. The sky was taken to be adjacent and to the south for RDl, and adjacent 
and to the north for BD3. As a check on this small beam photometry, aperture photometry 
was obtained of the combined flux from both objects, and this total was distributed between 
the two objects based on the one dimensional distribution of flux along the line between the 
two sources. The two methods agreed well, and the magnitudes are reported in Table 1. 
For RDl we measure J= 24.3 ± 0.3 mag and K> 23.1 mag (3(7 limit). For BD3 we measure 
J= 24.1 ± 0.3 mag and K= 23.0 ± 0.4 mag. In Table 1 we have also included the 1-band 
magnitudes of these objects as reported by Dey et al. (1998). 

To use the broadband magnitudes to derive continuum fluxes requires correction for 
contamination by strong emission lines contained within the fllter bandpasses. In the 
case of RDl the J fllter contains redshifted GUI] 1909A (at 1.206/xm), while the K fllter 
contains redshifted [OII]3727A (at 2.355/im). Heckman et al (1998) have recently compiled 
a composite UV spectrum of nearby starburst galaxies that shows the GUI] emission-line 
equivalent width to be generally less than lOA. This is small compared to the uncertainty 
in the measured J-band magnitude of RDl. Similarly, the average radio galaxy spectrum 
in McCarthy (1993) can be used as an upper limit on the strength of GUI] in RDl. If the 
equivalent width of GUI] in RDl is the same as in a typical radio galaxy (32A), the GUI] 
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line would contribute only about 8% of the total flux in the J filter. Thus we conclude that 
the correction for line emission in the J-band for RDl is not significant. However, unlike 
GUI], the [Oil] 3727A line can be quite strong in star forming galaxies, having rest-frame 
equivalent widths of 50 — 100 A (e.g. Cowie et al. 1995, Gallagher et al. 1989). At an 
equivalent width of lOOA [Oil] would contribute 14% of the (3a) limit we measure for the 
K-band flux from RDl. 

In the case of BD3, there are no strong hues expected to be present in the J filter, 
unless Mgll at 2800A has a broad, blue wing, which is unlikely for a starburst galaxy. The 
only emission feature that could contribute to the K-band measurement is the H/3 line. 
While potentially strong, the redshifted wavelength of H/? (2.44//m) places it at the very 
edge of the K filter and it should have a neghgible effect on the total flux. 



4. Discussion 

The most important result of these observations is that RD 1 was detected at J. This 
immediately implies that it is a much redder system than expected for any dust free, young 
galaxy model. The expected color of an unreddened, star forming galaxy with fi, ~const is 
I— J~ 0.5 mag and I— K~ 1.5 mag, while the apparent color of RDl is I— J= 2.2 ± 0.3 mag 
and I— K < 3.5 mag. Note, that source BD3 is also redder than expected for an unreddened 
starforming galaxy, having I— J= 1.2 ± 0.3 mag and I— K= 2.2 ± 0.5 mag. 

For comparison, the colors of the galaxy GL1358-I-62G1 at a redshift of ^ = 4.92 
(Pranx et al. 1997, Soifer et al. 1998) are I-J= 1.2 ± 0.1 mag and I-K= 2.1 ± 0.2 mag. 
For GL1358-I-62G1, Soifer et al. have shown that galaxy models with reddenings of 0.1 < 
E(B— V) < 0.4 mag provide substantially better fits to the overall energy distribution than 
do models with no reddening. 
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Similarly, Sawicki & Yee (1998) find that the rest- frame UV-optical colors of 
spectroscopically confirmed z > 2 Lyman break galaxies in the HDF imply significant 
internal dust obscuration at these redshifts. Again, young populations with E(B— V) ~ 0.3 
provide better fits to the data than do older, un-reddened models. Although RDl and BD3 
appear redder than the six, z > 3 sources in Sawicki & Yee, there is considerable scatter 
among the HDF galaxies. The range in I— J color of the HDF galaxies (closest to the J— K 
colors of RDl and BD3 in rest wavelength given the smaller redshifts of the Sawicki & Yee 
galaxies) is —0.98 <I— J< +0.65 mag, while the J— K colors of RDl and BD3 are < +1.2 
mag and +1.1 mag, respectively. Since we have argued above that the J-band magnitude of 
RDl is not likely to be significantly contaminated by strong line emission, it seems that the 
measured continuum colors of RDl and BD3 are likely to be affected by dust. 

With only an I— J color, and a limit on the I— K color, it is difficult to disentangle 
the effects of intrinsic galaxy color and reddening for RDl. However, we can identify a 
range of allowable ages and reddening values for a set of simple galaxy models by using the 
synthetic spectral energy distributions of Bruzual & Chariot (1993, 1996). Because we are 
simply trying to set the constraints on reddening and stellar age, we choose to fit aging, 
instantaneous starburst models, since these are the intrinsically reddest galaxies at a given 
age, and thus they require the least amount of dust reddening to match a given observed 
spectrum. Continuous star formation models will always require more reddening by dust. 

Before comparing the data to the models, a correction must be made to the observed 
I-band flux density for the presence of foreground absorption. Dey et al. (1998) correct 
their measured I-band magnitude for Lya emission and estimate 1= 26.5 + 0.1 mag in a 1.5" 
diameter beam (see Table 1). However there is no light detected blueward of Lya in the 
spectrum, and much of this absorption may occur at redshifts well below z = 5.3. Thus, the 
intrinsic rest-frame UV flux density of RDl is underestimated if the Lya emission- line alone 



-9- 



is taken into account. By comparing the Keck I-band filter and CCD responses to the RDl 
spectrum, we have estimated that this correction amounts to approximately 34% of the 
measured I-band light, assuming all of the light blueward of Lya is absorbed in intervening 
material. Once this is taken into account, RDl has an I-band magnitude of ~ 26.2 mag. It 
is this value that we use in all subsequent modelling. 

The comparison of the present data with the Bruzual & Chariot models shows that for 
RDl, either young, dusty models, or "old", dust-free models provide acceptable fits to the 
I, J, and K-band data. We have used both solar and 0.2Zq metallicity models, although 
we consider the latter more representative of the galaxies at these redshifts. Both sets 
of models require significant absorption. Ay > 0.5 mag, for stellar populations with ages 
significantly less than 10^ yrs. At an age of 10® yrs, the 0.2Zq model requires Ay = 0.3 
mag, while the solar metallicity model requires no reddening to fit the data. Both sets of 
models deliver unacceptable fits to the data for ages of 5 x 10® yrs or greater, regardless of 
the reddening. In all cases, the SMC reddening curve of Gordon & Clayton (1998) and the 
Bruzual & Chariot models with a Salpeter IMF have been used. 

If RDl is indeed a reddened, young galaxy then both its luminosity and star formation 
rate are larger than those estimated by Dey et al.(1998). If RDl has an age of significantly 
less than 10® yr and a metal abundance ~ O.2Z0, the visual extinction is larger than 
~ 0.5 mag, and the extinction at 1270A is a minimum of ~ 2.7 magnitudes, based on the 
reddening curve for the SMC (Gordon and Clayton 1998). Assuming the reddening to 
the Lya line is the same as that to the far UV continuum, the Lya and UV continuum 
luminosity are larger by factors of 10 — 15 than those determined by Dey et al. (1998) for no 
reddening. The monochromatic luminosity of RDl at 1500A is then a factor of ~ 3 larger 
than L^5QQ determined by Dickinson (1998) for z ^ 3 starburst galaxies. The star formation 
rate in RDl is then ~ 50 — lOOh^Q MQyT~^ , which is comparable to the extinction-corrected 
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values estimated by Sawicki & Yee for the z>2 Lyman break galaxies in the HDF. 
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Table 1: Photometry of RDl and BD3 



Object 


I J 


K 




mag mag 


mag 


RDl 


26.5±0.1^'''24.3±0.3 


23.r 


BD3 


25.1±0.P 24.1±0.3 


23.0±0.4 



"I-band magnitudes are from Dey et al. (1998) as measured in a 1.5" diameter circular beam. 

''The I-band magnitude of RDl has been corrected for Lya emission in the bandpass as per 
Dey et al. 

3 sigma limit, corresponding to a 1 sigma limit of 24.3 mag 
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Figure Captions 



Fig. 1.— Images of the field of 6C0140+32 RDl at J and K obtained witli tlie W.M. 
Keck Telescope. The scale and orientation of each field is shown. The images show the 
locations of the objects RDl (z — 5.34) and BD3 (z — 4.02) identified by Dey et al. 
(1998). After coadding the data in each filter, the images have been smoothed with 
a gaussian matching the average FWHM of a stellar source in the final mosaic (see 
text for details). The positions of RDl and BD3 are 10.2" west, 5.6" north and 10.6" 
west, 6.5" north of star A {a =01:43:43.67, 5 =+32:53:54.3 J2000), respectively, as 
given in Dey et al. 
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